Magneto-transport in impurity-doped few-layer graphene spin valve 
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Using Keldysh nonequilibrium Green's function method we study the spin-dependent transport 
through impurity-doped few layer graphene sandwiched between two magnetic leads with an ar- 
bitrary mutual orientations of the magnetizations. We find for parallel electrodes magnetizations 
that the differential conductance possesses two resonant peaks as the applied bias increases. These 
peaks are traced back to a buildup of a magnetic moment on the impurity due to the electrodes spin 
polarization. For a large mutual angle of the electrodes magnetization directions, the two resonant 
peaks approach each others and merge into a single peak for antiparallel orientation of the electrodes 
magnetizations. We point out that the tunneling magnetoresistance (TMR) may change sign for 
relatively small changes in the values of the polarization parameters. Furthermore, we inspect the 
behaviour of the differential conductance and TMR upon varying the temperature. 

PACS numbers: 85.75.-d,75.47.-m,71.55.-i 



I. INTRODUCTION 

In the past few years much effort was devoted to the 
investigations, fabrications and utilizations of graphene 
samples^— culminating in a series of fascinating find- 
ings such as the anomalous quantized Hall effect, the 
absence of weak localization and the existence of min- 
imal conductivity^. These phenomena underline the re- 
markable potential of graphene in future advances in 
nanoscience. Fueled by this development novel graphene- 
based technological applications and devices are cur- 
rently envisaged. In this respect, a special attention has 
been paid to the spin dependent transport in graphene 
and graphene-ferromagnet heterostructures, as prototyp- 
ical spintronic device^—. For example, Hill et alJ? fabri- 
cated graphene spin valves and observed a 10% change in 
the resistance as the electrodes magnetizations switch ori- 
entation from a parallel to an antiparallel configuration. 
Recent experiments on the spin injection in a single layer 
graphene show a rather long spin-flip relaxation length 
w 1/im at room temperature^. The spin injection into 
a graphene thin film has been successfully demonstrated 
by using nonlocal magnetoresistance measurements^—. 
Wang et al? measured the magnetoresistance of meso- 
scopic graphite spin- valve devices and observed a cusplike 
feature of the magnetoresistance versus the applied bias. 
Ding et a/.— studied theoretically the spin-dependent 
transport through the graphene spin valve device, and 
pointed out that a pronounced cusplike feature at zero 
bias is due to the result of a subtle combined effect of 
graphene and the conventional spin-valve properties. 

Recently, it has been demonstrated^ that adatoms can 
be precisely positioned on graphene. Numerous works 
evidence that the adatoms may create new many-body 
states in graphene and lead to extraordinary properties 
such as magnetism^— and Kondo effect* 2 ^ - — , which are 
different from the case of impurities in an ordinary metal. 
Therefore, the modification of the properties of graphene 
by the impurity atoms may also influence the magneto- 



transport in graphene nanojunctions, an issue which is 
addressed here. In this work, we study theoretically the 
spin dependent transport through few layers of graphene 
in the presence of impurity atoms. The method is based 
on the standard Keldysh nonequilibrium Green's func- 
tion approach, as described i n 26 i 27 . We find that if the 
electrodes spin polarization vectors are parallel, the ap- 
plied bias dependence of the differential conductance ex- 
hibits two resonant peaks signaling the formation of the 
impurity magnetic moment due to the electrodes mag- 
netization. With increasing the mutual angle 9 of the 
electrodes spin polarization vectors the two peaks develop 
gradually into a single peak at 6 = w. Thus the electrodes 
magnetization orientations may be used to switch on and 
off the magnetism of the impurity atoms. We also inves- 
tigate the dependence of the tunnel-magnetoresistance 
on 0, on the temperature, and on the spin polarization 
degrees of the electrodes in details. 

II. SINGLE LAYER GARAPHENE 




FIG. 1. A schematic illustration of the system considered in 
this work. A nanotunnel junction made of one or few layers of 
graphene with impurity atoms is connected to two magnetic 
leads. The magnetic moments of the leads are aligned by a 
relative angle 8, and the coupling matrix between a(a = L, R) 
electrode and graphene is T ka . 
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We consider a monolayer graphene sandwiched be- 
tween two ferromagnetic electrodes. Some impurity 
atoms are absorbed on the top of carbon atoms in the 
graphene sheet, as shown in FigfTJ The moment Mi of 
the left electrode is assumed to align along the y direc- 
tion, while the moment of the right electrode devi- 
ates from the y direction by a relative angle 9. A bias 
voltage V is applied between the left and the right elec- 
trodes. The electric current flows in the x direction. The 
Hamiltonian of this system reads 

H = H L + H R + H G + H T + H i + H f . (1) 

Here, H L (H R ) describes the left(right) electrode: 

h l = E e kL<r4 io .c kL(T (2) 

Hr = y^[gfl(k)-gMfl cos 0]<4 Ro .Ckfl <r -M,R sin flcj^Ckflq 

k,cr 

(3) 

where £kaer is the single electron energy, c ] { CQ0 .(ckacr) is 
the usual creation (annihilation) operator for an electron 
with the momentum k and the spin a in the a = L, R 
electrode. 

The tight-binding Hamiltonian of the electrons in 
graphene is given by 

H G = -t ]T (atAa+H.c), (4) 

where a} a {o,i,a) creates (annihilates) an electron with 
the spin a on the position of the sublattice A, 
b\ a .(bi j(T ) creates (annihilates) an electron with the spin 
a on the position R s ; of the sublattice B, t is the nearest 
neighbor((i, j}) hopping energy. In the momentum space 
Eq.Q can be rewritten as 

H G = ^[<Xq)a q > q - + <^(q)*6+ ff a qCT ] (5) 

q,cr 
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where </>(q) = -t E ^' Si with 5 i = § (h VS, 0), S 2 = 

i=l 

|(1, — -\/3, 0), S3 = o(l, 0, 0)(here a is the lattice spac- 
ing). The Hamiltonian ([5]) can be easily diagonalized 
and one can show E±(q) = ±i|</>(q)|, which can be lin- 
earized around the K points of the Brillouin zone and 
have a dispersion given by 

E ± (q) = ±v F \q\, (6) 

where vf = 3fo/2 is the Fermi velocity of an electron. 
The coupling between the electrodes and graphene is 
modeled by 

H T = 7= E [WLraq* + H - C '] . ( 7 ) 
kqacr 



where Tk Qq is the coupling matrix between the a elec- 
trode and the graphene; N is the number of sites on 
the sublattice A. Assuming some impurity atoms to be 
placed on the top of the sublattice A, for the simplicity, 
we neglect the correlated interaction between the impu- 
rities, and express the hybridization with the localized 
impurity states as 

Hi = 7^ E + a{M, (8) 

where f} a {fia) creates(annihilates) an electron with spin 
a at the impurity at the site I, and Vi is the random 
and strength distributions of the hybridization satisfying 
(ViVi')dis = Vq6u>, where (■■■)dis denotes the impurity 
averag e 28 ' 29 . 

Hf describes the impurities: 

H } =Y,^flji. + Un nnii , (9) 

Icr 

where n\ a = j\ a jia is the occupation number operator, 
£0 is the single electron energy at the impurity, and the 
Coulomb interaction is included by a finite U Ander- 
son term. For simplicity, we adopt a mean field ap- 
proximation to the electronic correlations at the impu- 
rity, Uniffiii ~ UJ2( n iv)fLfiv ~ u ( n n)( n il)i meaning 

a 

that the present theory is reliable for temperatures above 
the Kondo temperatures. The impurity Hamiltonian is 
rewritten as H f = £iaf} a fi<j with e i(T = e + U {nw). 

The electric current can be calculated from the time 
evolution of the occupation number operator of the left 
electrode. 

1 P 

I = e{tf L ) = ~{[H,Ar L ]), (10) 
where Ml = J2 c t.La CkLa '- Using the nonequilibrium 

kcr 

Green's function method, Eq. ([T0|) can be further ex- 
pressed as 

1 = - l iS^Tr{[Ql{e)-g-{e)]f L {e) + g<{e)}T L {e), 

(11) 

where TV is the trace in the spin space, / Q (e) is Fermi 
distribution function, Q r a {e) = £ (G qaq , a (e))d ls and 

qq' 

Ga ( £ ) = E( G q a,q'a( £ ))*s are 2 x 2 matrices representing 
qq' 

the retarded green's function and the lesser Green's func- 
tion respectively. In the calculation of Eq. pTj) . we as- 
sume that the dominant contributions to tunneling stem 
from the electrons near Fermi level, and hence assume 
the linewidth function to be independent of q. Thus, we 
have 
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The lesser Green function Q< (e) can be calculated by 
the Keldysh equation Q<(e) = g r a (e)^<(s)g^(e). To ob- 
tain T,<(e), we invoke Ng's ansatz^: E<(e) = T,^(e)B, 
where Ejf(e) = i[T L f L {e) + RT R I#f R (e)] is the lesser 
self-energy of the clean graphene system with 



R 



COS : 



sin ■ 



B is determined by the condition B = (Eg(e) — 
SgCe))- 1 ^- 1 ^)-^ 1 ^)) with ES( £ )(Eg(e)) denoting 
the retarded(advanced) self-energy of the clean graphene 
system. Under these considerations, we finally obtain 

1=1 J ^TrF L Ql{e)^9l(e)\\J R {s)- f L (s)}, (13) 

where E = RT R R^B. The remaining task is to calcu- 
late the retarded Green's function g r a (e). 

By the equation of motion and then impurity average, 
we can derive 



g r a (s)=g r a (e)[l-9 r a m r (e)}-\ 



(14) 



where gl< a (e) = ^ E Sqa,qa( £ ) and E r = n^g c - 
q 

|[rj,(e) + RY R (e)R)\ with the impurity density, 

5qa,qa( £ ) = ( £ ±i v )^- | 0(q) |2 ' aI1 d 



.9c 



e-e t 





Introducing a cutoff fe c leads to 

9a(e) = -Fo(e) ~ iirp (s), 



(15) 



(16) 



E | £ 2_ g 2| 

F o(£)-^2ln ^ 



Po (e) = ^9(D-\s\) (17) 



with £) = vpk c denoting a high-energy cutoff of the 
graphene bandwidth. k c is chosen as to guarantee the 
conservation of the total number of states in the Brillouin 
zone after the linearization of the spectrum around the K 
point, this is achieved following the Debye's prescription. 
We know from the selfenergy in Eq. (fT4| that introducing 
the impurities is behaving like changing the Fermi level 
with the impurity concentration. A similar effect, i.e. the 
Fermi level can be tuned from lying in the conduction 
band to the gap, is experimentally demonstrated by dop- 
ing Ca in a topological insulator Bi2Se3 31 . The position 
of the impurity level is also important since a resonance is 
clearly seen also as shown below. In Eq. (fT"3)) . we further 
set the symmetrical voltage division: (Xl,r — Ep ± \eV , 
and put Ep = in the numerical calculations. The oc- 
cupation of an impurity (n a ) which appears in the spin- 
dependent energy of impurities should be treated in a 
self-consistent way via the relation 



def(e)\mG r f (e), 



(18) 



where 

G r f (e) = (l-g^ig^-TairiiVi+g- 1 ^)}- 



(19) 



with E r d (e) = -|[r L (e) + RT R (e)R?]. In terms of the 
result of Eq. (fT"3]) . the TMR can be obtained according to 
the conventional definition 



TMR 



no) - i_M 

1(0) 



(20) 



where 1(0, it) is the current flowing through the system 
in the parallel (antiparallel) configuration. 



III. BI- AND TRI-LAYER GRAPHENE 

Let us consider a bilayer (trilyer) graphene by adding 
a second (third) layer according to Bernal-type stacking 
order (ABA). The inter-layer hopping energy is denoted 
by t p . By a straightforward generalization of the mono- 
layer case we obtain the retarded Green's function for the 
bilayer graphene as 



g r 2a (e)=A 2 (e)[l-A 2 (e)i: r (e)}-\ 



(21) 



where 



A 2 (e) 



JV S •uJ.|q|4-2£2 t ,J,| q |2 +e 4_ t 2 £ 2 ■ (22) 



For trilayer graphene we find 

G r 3a (e)=A 3 (e)[l~A 3 (e)J: r (e)}-\ (23) 

where 

Me) = -TrE- 
q 

with 



« |q| 6 +S2«||q| 4 + C 2 ^|q| 2 +D2 I 24 ) 



At = s, B 1 = t p e - 2e J 



d=e b - 2t z p e'\ B 2 = -3e z , C 2 = -2t z p e z + 3e 4 , 



and 



2_4 



Substituting Eqs.(|2Tj) and (|23J in Ea.(|T3j). we obtain the 
electric current through the bi- and tri-layer graphene 
with impurity atoms, and also the TMR according to 
Eq.(|2U|). The summation over q in Eqs.(|22]) and (|24|) 
may be performed by taking the continuum limit with a 
cutoff D and expanding the integrand in terms of partial 
fractions. The final results have the explicit form 



A 2 (e) 



1 f sgn(e)t 

D 2 t 2 



1 I (D 2 



D- 



W i. 2 — i (-U^-x 2 )x 1 
+ fln |(^x 1 )(^-x £ ) |} _, 

x[sgn(^L)6»(0 < x x < D 2 ) 
-sgn(^)0(O <x 2 < D 2 )} 
+ l[ sg n(^)9(0 < Xl <D 2 ) 
+sgn(^)0(O <x 2 < D 2 )]}, 



f sgn(e)t p 
I 2 



(25) 
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where 9(x) is the step function, and x-yp = e 2 ±i p |e|. For where 
A = (2S| - 9B 2 C 2 + 27 D 2 ) 2 + 4(-S| + 3C* 2 ) 3 > 0, 



A 3 (e) 



where 



1 

B 3 



r{[^i + 



j 4i(x 2 +3:3)3:i+-Bi.ti-Ai3:2^3+Ci 1 
(x 2 -xi)(x 3 -xi) J 



x (arctan 



arctan ■ 



D' 



\J X 2 X 3 -{X2 +X3 ) 2 /4 

-(x 2 +X3)/2 



•\/x 2 X3-(x 2 +X3) 2 /4 



(x 2 +x 3 )/2 



r) 



V 2:22:3 -(2:2 +a=3 ) 2 /4' 
Ai(x 2 +j:3)2:i+-Bi2;i-Ai2;2 3:3+Ci 



In 



(£C 2 — Xl)(x 3 — Kl) 

-D 2 — (2:2+2:3) +2:2 2:3 



2:22:3 
(x 2 +x 3 )/2-xi 



2^/x 2 x 3 - 

x (arctan 
+ arctan 

' dec 



(x 2 +x 3 ) 2 /4 

_D 2 -(x 2 +x 3 )/2 
\J X2X3 — (x 2 +x 3 ) 2 /4 

(2:2+2:3)72 

•y/ x 2 x 3 -(x 2 +2;3) 2 /4 



)]} 



^sgn(^)0(O < x x < # 2 )^f 



Ai (2:2 +2:3)2:1+512:1 — A 1 x 2 x 3 +Ci 
(x 2 — xi)(x 3 — xi) 



(26) 



x x = 

X2 = 

%3 = 

T = 



3 

B 2 
3 

B2 _ 
3 , 
2(B 2 



2V Bg-3C a cos( arccosT ); 

2\Z-B|-3C 2 , 2ir+arccosT \ 
3 Lub l 3 J: 

gf-3C 2 /47r+a£ccosT\ 
3 LU! H 3 J) 

-3C 2 )B2-3(B2C2-9-D2) 



(30) 



2(Bf-3C 2 )? 



rtcr), we need to calculate the im- 

i=2,3;/(* - = -**(* - 



To determine 

purity Green's function G crcr 

tf)({fa(t), fl, (t 1 )}) for the bilayer and trilayer graphene. 
Using the equation of motion method we find 

GVjie) = (l-A i ( £ )E5)[ 9c - 1 -A i ( £ )(n^ 2 + ffc - 1 S-)]- 1 

(31) 

Substituting Eq.([3"Tj) in Eq.(fI5|). we find (n a ) self- 
consistently. 



+ imH~ 2B 2 + 9^2C 2 - 27J 2 
V (2B 3 - 9B 2 C 2 + 27L> 2 ) 2 + 4(-B 2 + 3C 2 ) 3 }^ 

- ^{-2^1 + 9^2-27^2 

-y/(2B* - 9B 2 C 2 + 27D 2 f + 4(-5 2 + 3C 2 ) 3 }i 

(27) 



£2,3 = 



B2 
3 



-2S| 



95 2 C 2 - 27L> 2 



±\/(25 3 - 9B 2 C 2 + 

2i/3 1 zi3 2 



+*" 



27D 2 ) 2 +4(-B 2 + 3C 2 ) 3 }; 
+ 9B 2 C 2 - 27D 2 



TV( 2B 2 - $B 2 C 2 
For A < 0, 



27^ 2 



4(- 



SI 



3C 2 ) 3 }J 
(28) 



Aa(e) 



x In |^!=2 



[ii(i2+J:3)+flll^l- A 1 X 2 X3+Cl - 



(x 2 — Xl)(x 3 - 

fA 1 (x 2 +x 3 ')+B 1 



xi) 



[Ai(x 2 +x 3 )+_Bi]xi 



2:3-2:2 
-Aix 2 x 3 +Ci 



+ [ 



(x 2 — xi)(x3-x 2 ) J 
Ai(x2+x 3 )+Bi 
1 2:3—2:2 
IA 1 (x2+x 3 )+B 1 ]x 1 —A 1 x 2 X3+Ci 

-x 2 ) 

A1X2X3+C1 1 



In 



^1 



(x 2 — Xl)(x 3 - 

[Ai(x 2 +x 3 )+-Bi]xi- 



X In I 



(x 2 - 

£> 2 -x 3 



1} 



2:1) (2:3 -2:1) 
ttV? 



3 1 J D~ l ^ \ auj ' v 

[Ai(2: 2 +2;3) + Bi]2;i-Ai2;2 2;3+Ci 



f{sgn(fM#(0 < X! < D 2 ) 



(2:2 -Xl) (x 3 -Xl) 

_ sgn (^ )g( 0<, 2< ^2 )[ Ai(x 2 + ! 3)+Bi 

I [Ai(x 2 +x 3 )+i?i]xi— Aix 2 x 3 +Ci I 
' (x 2 -xi)(x 3 -x 2 ) 

+ S gn(tfM0<x 3 <^ 2 )[ Al( ^ +Bl 

I [Ai(x 2 +x 3 )+_Bi]xi-Aix 2 x 3 +Ci 
(x 2 -xi)(x 3 -x 2 ) 
[Ai(x 2 +x 3 )+_Bi]xi— Aix 2 x 3 +Ci 1 



(x 2 — xi)(x 3 -xi) 



-]}> 



(29) 



0.04- 




FIG. 2. (color online) The bias dependence of the differential 
conductance G for different polarizations P at T = O.OOiro 
(a) and for different temperatures T at P = 0.4 (b) for a single 
layer graphene junction and for parallel configuration of the 
electrodes magnetizations. The other parameters are taken 
as V f = 3r ,ni = 0.1, eo = 0.1F ,U = T ,D = 5F , where 
To stands for the coupling between the scattering region and 
the electrodes. 



5 



IV. NUMERICAL ANALYSIS 



To illustrate the nature of the spin-dependent trans- 
port we performed numerical calculations assuming the 
linewidth function (e) to be independent of the energy 
within the wide band approximation. Introducing the 
degree of the spin polarizations of the left and the right 
electrodes, Pl and Pr, and assuming that the two elec- 
trodes are made of the same material, i.e. Pl = Pr = P, 
we can write = r^j" = r (l ± P) where r describes 
the coupling between the graphene and the electrodes in 
absence of an internal magnetization, and is taken as the 
energy scale in the following numerical calculation. 

FigEla) shows the bias dependence of the differential 
conductance G = dl/dV in a single layer graphene sys- 
tem for the different P under the parallel configuration 
of the electrodes magnetizations. When P — (nonmag- 
netic electrodes), the differential conductance possesses a 
single resonant peak as the bias voltage increases (or de- 
creases) from zero up to few Tq (or — Tq). This peak cor- 
responds to a resonant tunneling through the renormal- 
ized spin-dependent impurity level Sq + Un a . Thereby, 
this situation reflects the fact that the impurity is not 
magnetized in the absence of the internal magnetization 
of the electrodes, i.e. = n^. When P^O, the afore- 
mentioned single peak in the conductance splits into two 
peaks indicating the positions of the resonant tunneling 
at £o + Un^ and Eq + E/nj, and signaling the form of the 
localized magnetic moment of the impurity. This be- 
havior is due to the fact that the exchange splitting of 
the density of state (DOS) in the ferromagnetic leads 
acts as an effective magnetic field with values well larger 
than those of an externally applied field 32 ' 33 . This ef- 
fective field acts on the graphene including the impurity 
atoms and influences strongly the spin states of the im- 
purities. With increasing P the peaks split further way 
from each other signifying that the polarization in the 
ferromagnetic leads enhances the localized magnetic mo- 
ment. With increasing P, the peak at the higher voltage 
is enhanced and that at the lower voltage is suppressed. 
This non-monotonic dependence of the peak amplitude 
on the polarization of the electrodes is completely differ- 
ent from the ferromagnet-quantum dot-ferromagnet(FM- 
QD-FM) system 3 -! - — ; its origin stems from the linear 
DOS of graphene. The dependence of the differential 
conductance on the bias for different temperatures in a 
parallel configuration is shown in Fig. [2fb) for a single 
layer graphene device. With increasing the temperature, 
the resonant peaks in the differential conductance de- 
crease and almost vanish at larger T. This temperature 
dependence of the conductance peaks is similar to that of 
a noninteracting, single-particle resonance in the multi- 
channel model. The mechanism for this destruction of 
the peak is that at high temperatures not all electron 
states of the low-lying subbands are fully occupied due 
to the occupation of the next subbands^. However, near 
V = 0, the conductance increases when the tempera- 
ture is raised. This characteristic feature is the same as 
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FIG. 3. (color online)The bias dependence of the differential 
conductance G for different energies eo of the impurity level 
at U = To, 9 = and rii =0.1 (a), for different interaction 
strengths U at eo = O.iro, 8 — and rii = 0.1 (b), and 
for different angles 9 at £o = O.lFo, U — To and n» = 0.1 
(c) as well as for the different impurity concentrations rii at 
£o ~ O.iro, U = To and 9 — (d). A single layer of graphene 
is considered. The other parameters are taken the same as 
those of FigH 



that in the clean undoped graphene systems^. The rea- 
son is the graphene acts as a barrier at the zero energy 
point. Near zero-bias voltage, the increase of the ther- 
mally excited electrons with the temperature enhances 
the conductance. 

The bias dependence of the differential conductance of 
a single layer graphene device for different Sq and differ- 
ent interaction U is shown in Fig. and (b). With 
increasing eo an d U, the resonant peaks in the conduc- 
tance shift towards large bias voltages, and the peaks am- 
plitudes become large which is also related to the linear 
DOS of the graphene. In particular, one can find that the 
interval between the conductance peaks increases with 
the interaction U, in contrast it decreases with Eq. This 
is because the DOS around the impurity energy level is 
enhanced as Sq grows, which decreases the localized mag- 
netic moment on the impurit y 16 ' 18 . FigJ3fc) shows for 
a single layer graphene tunnel junction the bias depen- 
dence of the differential conductance when varying the 
mutual angle 9 between the electrodes magnetization di- 
rection. With increasing 9, the split structure of the reso- 
nant peaks in the differential conductance is washed out, 
and eventually the two peaks merge into a single peak at 
9 = 7T, meaning that the impurity loses its magnetism in 
this situation. This behavior suggests that the magneti- 
zation electrodes can be used as a valve device to open 
or close the magnetism of impurity atoms in graphene. 
Similar phenomena has also been observed in the FM- 
QD-FM syste m 32 ' 38 , where the magnetic moment in the 
quantum dot is induced by different spin-dependent tun- 
neling rates between the two magnetic electrodes. In our 
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FIG. 4. (color online) The bias dependence of the TMR for 
different degrees of polarization P at T = O.OOiro (a) and 
for different temperatures T at P = 0.4 (b) for a single layer 
graphene tunnel junction. The other parameters are taken 
the same as those of Fig[2] 



case, the impurity atom immerses in the Dirac fermions 
sea. Its magnetic moment is turned on via a media- 
tion of the itinerant massless Dirac fermion not via a di- 
rect tunnelin g 16 ' 18 . This shows new features as discussed 
above and below. In Fig. GUd), we plot the differential 
conductance as an function of bias voltage for different 
impurity concentrations in a single layer graphene de- 
vice. One can easily observe that the resonant peaks di- 
minish with decreasing the impurity concentration, how- 
ever their positions do not change. This suggests that 
the magnetization is maintained although with descen- 
dent effect of the impurities on transport in graphene at 
lower concentration since the information of the magne- 
tizations of the electrodes is transmitted by the Dirac 
fermions in the same way. 

The bias dependence of the TMR, defined in Eq.®, 
for different polarizations P and the different tempera- 
tures T in a single layer graphene junction is shown in 
FigUJ A small hump at zero bias corresponds to an en- 
hanced magnetoresistance due to the result of a nontriv- 
ial combined effect of graphene and conventional spin- 
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FIG. 5. (color online) The bias dependence of the differential 
conductance G for different interlayer couplings t p at T = 
O.OOiro and P — 0.4 for the bilayer graphene system. The 
other parameters are taken the same as those of Fig(2] 



valve properties^. The existence of the impurity sup- 
presses the amplitude of this structure since the impu- 
rity influence on the DOS extends to the Dirac points. 
Additionally, it is found that the TMR as a function of 
the bias voltage exhibits two pronounced dips at lower 
and higher bias. The latter one is particularly interest- 
ing as the TMR there may change sign from positive to 
negative for small values of the polarization parameters. 
To account for this behavior one should take into account 
the angular dependence of the conductance from FigEJc) 
which shows clearly the change of the positions and the 
heights of the peak with the angle 9, leading to the ap- 
pearance of the negative TMR. For large bias (see Fig. 
EJb)), the graphene system with the impurity atoms be- 
haves like that in the absence of the impurities^, where 
TMR hardly changes with temperature T. This phe- 
nomena suggests that the presence of the impurity atoms 
does not affect the magnetoresistance at high bias volt- 
ages. However, for a low bias voltage, TMR displays 
many important features: With increasing temperatures 
the hump at zero bias diminishes since the increase in the 
electrical current from the contribution of thermally ex- 
cited electrons is faster for the antiparallel configuration. 
Additionally, the dip in TMR decreases with T. There- 
fore, due to their combinations, TMR near the zero bias 
voltage displays a broad peak instead of the "W" -shape 
feature with increasing the temperature T, and eventu- 
ally develops into a broad dip. 

Fig [5] shows the bias dependence of the differential con- 
ductance in a bilayer graphene system. When the in- 
terlayer coupling grows, the conductance at zero bias 
increases, while the resonant peak diminishes. The ex- 
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FIG. 6. (color online) The bias dependence of the differential 
conductance G for different layer graphene at T = O.OOiro, 
P = 0.4 and t p — 0.3Fq. The other parameters are taken the 



same as those of Fig[2] 



planation for this phenomenon is as follows: the quasi- 
particles in the bilayer graphene have a peculiar nature 
with features akin both to Dirac and to conventional 
fcrmions. The contribution of conventional fermions orig- 
inates from the interlayer coupling that supports a metal- 
lic bilayer graphene. Thus the increase of the interlayer 
coupling enhances the zero bias conductance in the bi- 
layer graphene. Furthermore, the metallic behavior due 
to the interlayer coupling enhances the screening of the 
impurit y 40 ' 41 , which leads to the decrease of the trans- 
mission probability of the electrons, thus suppressing the 
resonant tunneling through the energy level £o + n a U . 
The bias dependence of the differential conductance for 
the different layer graphene devices is shown in Fig. [5] 
The differential conductance for all the graphene layers 
displays two peaks corresponding to the resonant tunnel- 
ing through the energy level e + n a U . It is suggested 
that the impurity magnetic moment also exists for the bi- 



layer and trilayer graphene when the ferromagnetic elec- 
trodes are introduced. However, the distance of the two 
peaks decreases with increasing the number of graphene 
layers. This result stems from the fact that the impu- 
rity magnetic moment decreases with increasing the layer 
number—. 



V. SUMMARY 

In conclusion, using Keldysh's nonequilibrium Green's 
function method, we study theoretically the spin depen- 
dent transport through a tunnel junction made of few 
layers of graphene with impurity atoms. It is found that 
when the electrodes magnetization vectors are parallel, 
the bias dependence of the differential conductance ex- 
hibits two resonant peaks as the bias voltage increases. 
We assign this behavior to the formation of a magnetic 
moment on the impurity due to the electrodes spin polar- 
ization. With increasing the mutual angle 9 between the 
vectors of the two electrodes magnetization, the distance 
of the two resonant peaks in the differential conductance 
diminishes, and they merge eventually into a single peak 
at 9 = 7T. The reason is that the induced magnetic mo- 
ment on the impurity decreases with 9 and vanishes for 
9 = tv. This features may be exploited to develop a mag- 
netic valve device that operates based on the presence 
or the absence of the magnetism on the impurity atom. 
Furthermore, due to the ^-dependence of the impurity 
magnetic moment, the TMR may change sign for small 
values of the polarization parameters. 



ACKNOWLEDGMENTS 

The work of K.H.D. and Z.H.Z. was supported by the 
National Natural Science Foundation of China (Grant 
Nos. 10904007 and 60771059), the Natural Science Foun- 
dation of Hunan Province, China (Grant No. 08JJ4002 ), 
and the construct program of the key discipline in Chang- 
sha University of Science and Technology, China. J.B. 
and Z.G.Z. were supported by the cluster of excellence 
" Nanostructured Materials" of the state Saxony- Anhalt, 
and the DFG, Germany. 



1 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. 
I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. 
Firsov, Nature 438, 197 (2005). 

2 Y. Zhang, Y. W. Tan, H. L. Stormer, and P. Kim, Nature 
438, 201 (2005). 

3 A. K. Geim and K. S. Novoselov, Nature Materials 6, 183 
(2007). 

4 A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, A. K. Geim, Rev. Mod. Phys. 81, 109 (2009). 

5 E. W. Hill, A. K. Geim, K. Novoselov, F. Schedin, and P. 
Black, IEEE Trans. Magn. 42, 2694 (2006). 



N. Tombros, C. Jozsa, M. Popinciue, H. T. Jonkman, and 
B. J. Van Wees, Nature 448, 571 (2007). 

7 S. Cho, Yung-Fu Chen, and M. S. Fuhrer, Appl. Phys. 
Lett. 91, 123105 (2007). 

8 M. Ohishi, M. Shiraishi, R. Nouchi, T. Nozaki, T. Shinjo, 
and Y. Suzuki, Jpn. J. Appl. Phys. 46, L605 (2007). 

9 W. H. Wang, K. Pi, Y. Li, Y. F. Chiang, P. Wei, J. Shi, 
and R. K. Kawakami, Phys. Rev. B 77, 020402(R) (2008). 

10 K.H. Ding, Z.G. Zhu, and J. Berakdar, Phys. Rev. B 79, 
045405 (2009). 

11 S.Krompiewski, Phys. Rev. B 80, 075433 (2009). 



8 



12 K.H. Ding, Z.G. Zhu, and J. Berakdar Europhysics Letters 
88, 58001 (2009) 

13 J.C. Chen et al. , J. Phys.:Condens. Matter 22, 035301 
(2010) 

14 K.-H. Ding, G. Zhou, Z.-G. Zhu, J. Berakdar, J. Phys.: 
Condens. Matter 20, pp 345228 (2008). 

15 J.C. Meyer et al, Nature 454, 319 (2008). 

16 B. Uchoa, V. N. Kotov, N. M. R. Peres, and A. H. Castro 
Neto, Phys. Rev. Lett. 101, 026805 (2008). 

17 A. H. Castro Neto, Valeri N. Kotov, Johan Nilsson, Vitor 
M. Pereira, Nuno M. R. Peres, Bruno Uchoa, Sol. Stat. 
Comm. 149, 1094 (2009). 

18 K.-H. Ding, Z.-G. Zhu, J. Berakdar, J. Phys.: Condens. 
Matter 21, 182002 (2009). 

19 B. Uchoa, L. Yang, S.-W. Tsai, N. M. R. Peres, A. H. 
Castro Neto, Phys.Rev. Lett. 103, 206804 (2009). 

20 M. Hentschel and F. Guinea, Phys. Rev. B 76, 115407 
(2007). 

21 B. Dora and P. Thalmeier, Phys. Rev. B 76, 115435 (2007). 

22 P. S. Cornaglia, G. Usaj, and C. A. Balseiro, Phys. Rev. 
Lett. 102, 046801 (2009). 

23 H. B. Zhuang, Q. F. Sun, and X. C. Xie, Eur. Phys. Lett. 
86, 58004 (2009). 

24 Z. G. Zhu, K. H. Ding, and J. Berakdar, Europhys. Lett. 
90, 67001 (2010). 

25 T. O. Wehling, H. P. Dahal, A. I. Lichtenstein, M. 
I. Katsnelson, H. Manoharan, and A. V. Balatsky, 
larXiv:0906.2426l 

26 H. Haug and A. P. Jauho, quantum kinetics in transport 
and optics of semiconductors (Springer, Berlin, 1998). 

27 J. Rammer, quantum transport theory (Westview Press, 
2004). 



28 N. A. Sinitsyn, J. E. Hill, H. Min, J. Sinova, and A. H. 
MacDonald, Phys. Rev. Lett. 97, 106804 (2006) 

29 T. S. Nunner, N. A. Sinitsyn, Mario F. Borunda, V. K. 
Dugaev, A. A. Kovalev, Ar. Abanov, C. Timm, T. Jung- 
wirth, J. Inoue, A. H. MacDonald, and J. Sinova, Phys. 
Rev. B 76, 235312 (2007). 

30 T. K. Ng, Phys. Rev. Lett. 76, 487 (1996). 

31 J. G. Checkelsky, Y. S. Hor, M.-H. Liu, D.-X. Qu, R. J. 
Cava, and N. P. Ong, Phys. Rev. Lett. 103, 246601 (2009). 

32 A.N. Pasupathy, et al., Science 306, 86 (2004). 

33 Z. G. Zhu, Phys. Lett. A 372, 695 (2008). 

34 N. Sergueev, Q.-F. Sun, H. Guo, B. G. Wang, and J. Wang, 
Phys. Rev. B 65, 165303 (2002). 

35 H.-F. Mu, G. Su, and Q.-R. Zheng, Phys. Rev. B 73, 
054414 (2006). 

36 Z.-G. Zhu, G. Su, Q.-R. Zheng, and B. Jin, Phys. Rev. B 
70, 174403 (2004). 

37 B.J. van Wees, L.P. Kouwenhoven, E.M.M. Willems, 
C.J.P.M. Harmans, J.E. Mooij, H. van Houten, C.W.J. 
Beenakker, J.G. Williamson, C.T. Foxon, Phys. Rev. B 
43, 12431 (1991). 

38 P. Zhang, Q. K. Xue, and X. C. Xie, 
|arXiv: cond-mat /020 1465 1 

39 J. Nilsson and A. H. Castro Neto, Phys. Rev. Lett. 98, 
126801 (2007). 

40 D. P. DiVincenzo and E. J. Mele, Phys. Rev. B 29, 1685 
(1984). 

41 N. M. R. Peres, F. Guinea, and A. H. Castro Neto, Phys. 
Rev. B 73, 125411(2006). 



